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This report represents some new experimental results and the quantitative model describing large 
magneto-strain effect and main mechanical and magnetic properties observed in several ferromag- 
t^-j- . netic shape-memory alloys. The model application to giant magneto-strain effect recently found in 

■ some non-stoichiometric Ni-Mn-Ga alloys is discussed. 

o 

I. INTRODUCTION 

Some ferromagnetic shape memory alloys were recently suggested as a general way for the development of a new 
class of the magnetic-field-controlled actuator materials that will allow control of large strain effect by application of 
a magnetic field. Numerous candidate shape memory materials were explored including Ni 2 MnGa Co2MnGa, FePt 
CoNi, and FeNiCoTi during past few years. Large magnetically driven strain effects are expected have to occur in 
these systems. References [8-24] contain a detailed information on this subject. Up to date, the largest magneto-strain 
effects were achieved in Ni-Mn-Ga ferromagnetic shape memory alloys [1-7]. The best results on the field-induced 
strain reported short time ago [2] never exceeded 0.2-0.3% strain value in a magnetic field of order 8 kOe in martensitic 
state for nearly stoichiometric single-crystal samples of Ni 2 MnGa alloy. Recently, several research groups [6,7] have 
reported on the observation of super-large more than 5% magneto-strain effect in some non-stoichiometric Ni-Mn-Ga 
alloys close to that 5.78% value expected from the tetragonality aspect ratio of the martensite crystal lattice. New 
Ni-Mn-Ga alloys showing giant magneto-strain effect display simultaneously few interesting physical effects and new 
behavior for some magnetic and mechanical properties which is very different from that earlier observed in Ni2MnGa 
t^J- \ with lower (< 0.3 A) magneto-strain effect. 

■ One the new non-stoichiometric Ni4gMn3oGa22 alloy has recently been studied in details. This alloy undergoes the 
martensitic transformation at about 35C° from the L21 ordered Heusler type cubic (a = 0.582 nm) crystal structure 
into a tetragonally distorted structure with crystalline lattice parameters: a=b=0.594 nm and c=0.562 nm. The 
martensitic phase is ferromagnetic one and internally twinned. It contains a mixture of twin variants twinned usually 
on {110} planes. The orientation of their c-axes is nearly parallel to [100] directions of the parent cubic phase. 
Tetragonal symmetry direction is simultaneously a single easy magnetization axis for each martensitic variant. Very 
low 2-3 MPa compressive twinning stress along [100] directions needed to transform the multi-variant samples into a 
single-variant state has been found for Ni-Mn-Ga alloys displaying giant field-induced deformation effect. It is much 
less than it was observed in some earlier studied alloys (15-20 MPa [8,9]). As a result, new alloys which usually have 

O ■ very simple two-variant twinning microstructure and very low twinning stresses can be easily transformed between 
J-^ ' two single twin variants by application both the mechanical stress and the magnetic field through the mechanism of 
^ , twin boundaries motion. 

For instance, Fig.l shows an example of giant 5% magneto-strain effect associated with the direct and reverse 
magnetic- field-induced transformation between two single twin variants in Ni4sMn3oGa22 sample. The magneto- 
strain value achieved in this new alloy is insufficiently less compared to the natural crystallographic limit 1-c/a = 
5.78% expected for twinning in this material. 

The next Fig. 2 represents in situ observation of the transformation process between two twin variants of martensitic 
phase in the same Ni-Mn-Ga sample through the mechanism of the magnetic-field-controlled twin boundary motion. 
These experiments directly confirm the magnetic-field-induced twinning as the main physical mechanism responsible 
for a super-large magneto-strain effect found in this material. 

Simultaneously, some important steps were made during the past few years concerning the physical nature of the 
magneto-mechanical phenomena observed in ferromagnetic shape memory alloys. James and Wuttig [4] developed the 
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models based on a constrained magnetization orientation approach. O'Handley [5] first included the rotation magne- 
tization effects together with usual Zeeman energy term into consideration. Likhachev and Ullakko [3] have developed 
a general thermodynamic approach based on the Maxwell's relationships and derived the main magneto-mechanical 
equation representing both the magnetic and the mechanical driving force balance in ferromagnetic shape memory 
materials. The last model has been also successfully applied to explain quantitatively the results of experimental study 
of large magneto-strain effects in N^MnGa earlier obtained in [2]. In particular, according to our model estimation 
made in [3] a giant magneto-strain values can be achieved only in materials with very low (< 2 MPa) twinning stress 
values. Mechanical testing results obtained for giant magneto-strain Ni-Mn-Ga samples have completely confirmed 
these estimations. 

The present report represents some experimental results and several new effects recently found in some non- 
stoichiometric Ni-Mn-Ga alloys that are directly connected with their giant magneto-strain response. We also represent 
some new model background based on the magnetization free energy consideration that allows better understand its 
important role. Naturally, this approach gives the same final equations derived in our preceding publications. Finally, 
the model equations are applied to explain some new experimental results. 



II. MAGNETIZATION FREE ENERGY AND MAGNETIC DRIVING FORCES 

Perhaps, starting from a simplest physical idea first clearly formulated in [1], it is generally believed that a large 
macroscopic mechanical strain induced by the magnetic field in some ferromagnetic shape memory alloys is usually 
realized trough the twin boundaries motion and redistribution of different twin variant fractions in a magnetic field. 
In this case, the main thermodynamic driving forces should have a magnetic nature and be connected with high mag- 
netization anisotropy and significant differences in magnetization free energies for different twin variants of martensite 
[1-7]. It should be noted that some early models have considered each twin band as a single magnetic domain with 
a constant magnetization value and constrained magnetization direction. This approach predicted some unlimited 
linear growth of the magnetic driving force in a magnetic field applied. Respectively, it was expected that one can 
easily obtain significant magneto-strain effects in ferromagnetic shape memory alloys trough the twin boundary mo- 
tion simply increasing magnetic field value as high as needed. It has been understood later [5, 3] that it is not so, 
because each twin band is not a single magnetic domain but should contain the multiple internal magnetic domain 
microstructure and typical expected for the uniaxial ferromagnet 180° domain wall configurations. As a result, the 
local average magnetization inside of twins is not a constant and can be changed in a magnetic field due to both the 
magnetic domain wall motion and the rotation of magnetization. The most important conclusion following from this 
point of view is that the magnetic driving force can not be increased infinitely in a magnetic field but always limited 
its maximal value independent on the magnetic field and proportional the magnetic anisotropy constant. This value 
is low enough, so only the very soft (in sense of twinning stress) materials can be good candidates to show really high 
magnetic field induced strain effects. Therefore, the investigation of magnetic properties is the first important task 
for the magnetic driving force calculation and understanding of the magneto-mechanical behavior of Ni-Mn-Ga and 
similar type alloys. 

Some results showing main magnetic properties of Ni4§Mn3oGa22 in martensitic state are indicated in Figs. 3-4. 
First, the magnetization behavior of a single variant of martensitic phase has been studied. For this aim the internally 
twinned sample was initially transformed into a single variant state by application of compressive stress (> 3MPa) 
and then hardly constrained without unloading in a special holder. Such a procedure was needed to avoid the magnetic 
field induced motion of twin boundaries and resulting transformation between the different twin variants during the 
magnetization measurements. Two magnetization loops displaying field dependence of the magnetization m a (ft) 
along the easy (tetragonal symmetry) axis and also along the hard magnetization direction mt (ft) are shown in Fig. 3. 
Different mechanisms of magnetization for easy and hard magnetization directions are expected, respectively. In 
particular, 180° magnetic domain walls motion is expected to be the main magnetization mechanism for easy axis. So 
as, the rotation of magnetization without domain walls motion seems to be responsible for the magnetization behavior 
in a hard direction [5]. Corresponding magnetization free energies can be calculated from these magnetization data 
as functions of the magnetic field applied as follows: 

h h 
9a (ft) = - J m a (ft) dh; g t (ft) = - J m t (ft) dh; (1) 
o o 

These results are shown in Fig. 4. The energy of magnetic anisotropy defined as a difference of these free energies 
is also indicated on this plot. In particular, it can be easily found that the magnetic driving force responsible for 
twin boundaries motion is practically equal to the magnetic anisotropy energy. For this aim one should also know 
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the magnetization free energy in the multi-variant sate. Its dependence on the relative volume fractions of different 
martensite variants is especially important. 

Fig. 5, for instance, schematically shows a typical for some new Ni-Mn-Ga alloys two-variant twin microstructure, 
orientation of easy and hard magnetization axes for both twin variants and also the magnetic driving force alignment 
that usually give the maximal value of the magnetic field-induced strain in this material. In this case of simple two- 
variant twinning geometry the easy magnetization axis of the first variant (white area) is parallel to a magnetic field 
applied. For the second one (gray color) the magnetic field is applied in a transversal hard magnetization direction. 
Further, both these variants are called axial (a) and transversal (t), respectively. Taking the specific magnetization 
free energies for both twin variants g a (h) and gt (h) found, for instance, from the experiment one can write the total 
magnetization free energy per unit volume of the material as follows: 

g mag (h, x) = xg a (h) + (1 - x) g t (ft) (2) 

Therefore, due to magnetic anisotropy the magnetization free energy becomes dependent on the relative volume 
fractions x and 1 — x occupied by axial and transversal twin variants, correspondingly. According to our previous cal- 
culation of the magnetization free energies shown in Fig. 4 g a (h) < g t (ft). So, one can decrease from its maximal value 
at x — to the minimal one at x = 1 moving the twin boundaries as indicated in Fig. 5. The corresponding magnetic 
driving force moving twin boundaries along their normal directions can be found from the general thermodynamic 
rule as follows: 

= gt (ft) - g a (h) (3) 

h 

Therefore, a non-zero magnetic driving force responsible for twin boundary motion appears as a result of uniaxial 
magnetic anisotropy of Ni-Mn-Ga. This force is applied along the twin boundary normal direction and equal to a 
difference in magnetization free energies between the different twin variants. The magnetic driving force is dependent 
on the magnetic field applied and never can exceed its saturation level as shown in Fig. 4. The maximal estimated 
value about of 0.13 * 10 6 N/m 2 is achieved at ft > 0.8T for Ni4gMn3oGa22 alloy. Similar way, twin boundaries can be 
also driven by the mechanical stress applied. The corresponding normal mechanical driving force as a function of the 
applied stress a is known can be found as follows: 

fmec (<?) = S d — S ((J xx - CT yy) (4) 

Here, £o is the strain matrix associated with twinning transformation between two single variants of martensitic phase. 
In the XY reference system shown in Fig. 5 it is represented as a diagonal matrix as follows: 

f e \ 

eo = -e (5) 
\ 0/ 

where, e = 1 — c/a = 0.0579 is the twinning transformation strain value that can be easily estimated from the 
martensitic phase lattice parameters data: a ~ b — 0.594nm and c = 0.562nm. This is a maximal strain allowed by 
crystallography of twinning in Ni-Mn-Ga. 
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III. UNIVERSALITY RULES AND RELATIONSHIP BETWEEN THE STRESS-INDUCED AND 
MAGNETIC-FIELD-CONTROLLED TWINNING STRAINS 

No general way has been found up to date that would allow calculate the macroscopic twinning strain as a function 
of driving force. Nevertheless, one can expect that independently on the physical nature of these forces (magnetic, 
mechanical, etc.) the same driving forces applied to twin boundaries must produce the same macroscopic deformation 
effects. This evident physical requirement or universality rule denotes that the macroscopic twinning strain must be 
some universal function dependent only on the driving force value and independent on the physical source of this 
force. So, in case of plane twinning both non-zero components of the strain can be written through a single universal 
function e u (f) dependent on the normal driving force / as follows: 

e yy = Sxx = e u (/) (6) 

Using this equation one can represent both the mechanical and magnetic strain simply assuming / to equal to fmec (<?) 
or fmag (ft) respectively: 
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e ™«= (a) = -e™ c (a)=e u (e a) (7) 

and 

e™° (ft) = -e^ (ft) = e u (g t (ft) g a (h)) (8) 

It is not difficult to eliminate the unknown universal function e u (/) from Eqns.(7,8) and find that e' na9 (h) and 
e*™ ec ((j)will take equal values if only the corresponding mechanical and magnetic driving forces are also equal. In 
other words: 

£™ a9 (h) = e™^ (ct) , if : f mec (ct) = f mag (h) (9) 

For instance, in case of the uniaxial compressive stress (a xx — a zz = 0; <r yy = — ct) the mechanical and magnetic 
driving forces become equal at: 

o = <7 mag (ft) = £o 1 (g t (ft) - g a (ft)) (10) 
This immediately gives the following important relationship: 

e ma9 (h) = e mec (a mag (h)) (11) 

which allows performing the quantitative calculations of the magnctic-ficld-controllcd deformation effects in different 
ferromagnetic shape memory alloys by using the corresponding mechanical testing data. According to Eq.(ll) the 
deformation effect of a magnetic field is equivalent to some additional uniaxial compressive stress applied. This 
equivalent magnetic stress a mag (h) can be easily calculated from the Eq.(10). For this aim one needs only to know 
the magnetization properties of the material and use the magnetic driving force calculation data indicated in Fig.3- 
4. The maximal equivalent magnetic stress that can be developed in a magnetic field is at about 2.25 MPa in 
Ni4 8 Mn 3 oGa 2 2. 

On the other side, according to results of the mechanical testing experiments shown in Fig. 6 the uniaxial mechanical 
stress required to transform completely Ni4§Mn3oGa22 samples from one single variant of martensitic phase to another 
one through the stress-induced twin boundaries motion does not exceed 2.5 MPa. The maximal compression strain 
that can be achieved on this load is about 5.8%. Therefore, due to very low twinning stress value recently found in 
some Ni-Mn-Ga alloys the calculated magnetic stress is completely enough to achieve very high 5% magneto-strain 
value in these new materials. 



IV. MODEL CALCULATIONS 



Some calculation results that follows from our model consideration are presented in this section. In particular, 
Fig. 7 represents both the strain-stress hysteresis loop found from the experiment and also a set of model fitting 
curves corresponding to different values of twinning stress. Simple Fermi-like distribution functions were used as an 
appropriate fitting basis for analytic interpolation of the mechanical testing results. 

^ ec (-)= £ o(l + cxp(^-^))" (12) 

Here, ± denote loading and unloading curves respectively. cto and Act are characteristic stress parameters that can be 
associated with the start (<r s — a® — 2 Act) and finish (ct/ = cto + 2 Act) twinning stress values. The best fit between 
the model function and experimental data is achieved at cto = 1.56 MPa and Act = 0.26 MPa and indicated in Fig. 7 
by the solid circle labelled line. Additional mechanical hysteresis loops (a — b, c — e) shown here were also simulated 
to study the effect of different twinning stress values on the magneto-strain behavior. 

The corresponding model calculation results performed in accordance with Eq.(ll) and the experimental mea- 
surement data for the magnetic-field-controlled strain behavior in Ni4 8 Mn 3 oGa22 are plotted altogether in Fig. 8. As 
follows from these calculations the increase of twinning stress may cause the gradual decrease of the magnetic field 
induced strain from its maximal possible value 5.8% at ct/ < 1.67 MPa to 2% at ct/ = 3.32 MPa. In particular, 
one can observe the reasonable quantitative agreement between the 5.5% value expected from the model calcula- 
tion and 5.2% experimental one. Therefore, one can conclude that the low twinning stress value is a very critical 
physical parameter controlling a super-large magneto-strain effect in Ni-Mn-Ga. In particular, this can explain con- 
siderably lower magneto-strain effect found earlier in some nearly stoichiometric alloys which had at almost ten times 
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higher (15 — 20 M Pa) twinning stresses. So, the magnetic driving forces (that can produce no more than 2 — 3MPa 
equivalent stress) were only able insufficiently redistribute the relative twin variant fractions in these materials. As a 
result no more than 0.2% deformation effect was observed. 

The next Fig. 9 represents another important application of the model concerning one interesting physical effect that 
has been found in magnetization behavior. This effect is also directly connected with the magnetic-field-controllcd 
twin boundary motion found in Ni-Mn-Ga alloy. Some very special magnetization hysteresis loop is observed if the 
magnetic field is applied to a single variant of Ni-Mn-Ga alloy along its hard magnetization direction. At the initial 
stage the magnetization curve exactly coincides with a corresponding hard magnetization curve shown in Fig. (3). At 
higher magnetic field the sample is transformed by twinning until the second most appropriate twin variant with easy 
axis parallel to a magnetic field applied will be formed. During the next magnetic cycles the magnetization curve will 
always follow the easy magnetization path practically without any hysteresis. The quantitative calculation of this 
effect can be done according to the following magnetization equation: 



that has been first found in [3] and can also be easily obtained in the framework of our present approach by using the 
well known thermodynamic definition: 



and evident relationship x = e/eq between the strain and transformed fraction. Resulting behavior of the magne- 
tization can be found from the Eq.(13) and previously calculated field dependence of the magneto-strain, assuming 
there e = e mag (h). Calculations were also done for different twinning stress values. This effect shows also a good 
agreement between the model and experiment. Finally, one should make the following important conclusions: 



1. A super- large deformation effect at about 5% that can be induced by application of magnetic field less than 
0.8T in some new ferromagnetic non-stoichiometric Ni-Mn-Ga shape memory alloys seems to be the most interesting 
property of this material. 

2. Ni-Mn-Ga alloys displaying giant field-induced deformation effects have always very simple practically two- 
variant twin microstructure. As a result, these alloys can be easily deformed by application of the low < 2MPa 
mechanical stress or magnetic field less 0.8T through the mechanism of twin boundary motion. In particular, the 
reversible magneto-mechanical transformation cycle between two twin variants of martensitic phase can be induced 
by subsequent application of magnetic field along two different perpendicular directions. 

3. The magnetic driving force applied to twin boundaries is equal to a difference in magnetization free energies be- 
tween the different twin variants of martensite. This difference also characterizes the energy of uniaxial magnetization 
anisotropy and can be calculated from the corresponding magnetization measurements. The magnetic driving force 
achieves its maximal value 0.13 MN/m2 in a magnetic field higher 0.8T. In agreement with our earlier estimations 
and present model calculations this value is completely enough to explain high 5% magnetic field induced strain in 
some Ni-Mn-Ga alloys which have the very low < 2MPa twinning stress. 

4. There is a definite analogy between the deformation effects caused by the mechanical and magnetic driving 
forces. For instance, in both cases the macroscopic twining strain driven by the mechanical stress or the magnetic 
field applied can be expressed through the same universal function dependent on the corresponding (mechanical or 
magnetic) driving force. This universality rule allows to perform the quantitative calculations of the magnetic field 
induced strain as a function of field by using the mechanical testing results. 



[1] K. Ullakko J. Met. Eng. Perform., 5, (1996) 405. 

[2] K. Ullakko, J. K. Huang, C. Kantner, R. C. O'Handley and V. V. Kokorin, Appl. Phys. Lett., 9, (1996) 1966. 

[3] A. A. Likhachev and K. Ullakko, EPJdir, B2, (1999) pp. 1-9; Eur. Phys. J, B14,(2000) 263. 

[4] R. D. James and M. Wuttig, Philos. Mag. A 77, (1988) 1273. 



m (e, h) = m t (h) + (e/e ) {m a (h) - m t (h)) 



(13) 




(14) 



V. CONCLUSIONS 



5 



[5] R.C. O'Handley, J. Appl. Phys., 83, (1998) 3263. 

[6] K. Ullakko, A.Sozinov, and P.Yakovenko, E-print, ArXiv.org/cond-mat/0004211. 

[7] S.J. Murrey, R. Hayashi, M. Marioni, S.M. Allen, R.C. O'Handley., In Smart Structures and Materials 1999: Smart Materials 
Technologies, ed. Manfred Wuttig, USA. 1999, SPIE3675, pp. 204-211. 

[8] V. V. Kokorin and V. V. Martinov Fiz. Met. Metalloved., 9, (1991) 106. 

[9] V. V. Martinov and V. V. Kokorin J. Phys. (France) III 2, (1992) 739. 
[10] V. A. Chernenko, E. Cesari, V. V. Kokorin and I. N. Vitenko, Scr. Metall. Mater. 33, (1995) 1239. 
[11] P. J. Webster, K. R. A. Ziebeck, S. L. Town and M. S. Peak, Philos. Mag. B 49, (1994) 295. 
[12] R. Tickle and R. D. James, J. Magnet. Magn. Mat., (submitted) 

[13] K. Ullakko, P. T. Jakovenko, and V. G. Gavriljuk, New developments in actuator materials as reflected in magnetically 
controlled shape memory alloys and high-strengh shape memory steels, Proceeding Conf. on Smart Structures and Materials, 
ed. V.V. Varadan and J. Chandra, vol. 2715, pp. 42-50, SPIE, San Diego, USA, 1996. 

[14] R.D. James and M. Wuttig, Alternative smart materials, Proceeding Conf. on Smart Structures and Materials, ed. V.V. 
Varadan and J. Chandra, vol. 2715, pp. 420-426, SPIE, San Diego, USA, 1996. 

[15] James R.D. International Journal of Solids and Structures (UK), 37, (2000), 239. 

[16] R.W. Overholser, M. Wuttig, D.A. Neumann, Scripta Materialia, 40, (1999), 1095 

[17] K. Ullakko. Large-stroke and high-strength actuator materials for adaptive structure. Proc. Third International Conf. on 

Intelligent Materials 1996, ed. P. F. Gobin and J. Tatibouet, Lyon, France. 1996, SPIE2779, pp. 505-510. 
[18] K. Ullakko, J. K. Huang, V.V. Kokorin, R. C. O'Handley, Scr. Materialia, 36, (1997), 1133. 

[19] R. D. James, R. Tickle, M. Wuttig. Large field-induced strains in ferromagnetic shape memory materials. Int. Conf. on 

Martensitic Transformations ICOMAT-98, p. 1.14, San Carlos de Bariloche, Argentina, 1998. 
[20] R. Tickle, R. D. James, T. Shield, M. Wuttig, V. V. Kokorin. IEEE Transaction on Magnetics, 35, (1999) 4301. 
[21] V. V. Kokorin, V. V. Martinov, V. A. Chernenko. Scr. Met. et Mat., 26, (1992), 175. 
[22] V. V. Martinov. Journal de Physique IV France, Colloque C8, 5, (1995) 91. 

[23] V. A. Chernenko, C. Segui, E. Cezari, J. Pons, V. V. Kokorin. Journal de Physique IV France, Colloque C5, 7, (1997), 137. 
[24] S. Wirth, A. Leithe- Jasper, A. N. Vasil'ev, and J. M. D. Coey. J. of Magnetism and Magnetic Materials, , 167, (1997), 
L7-L11 



6 



Figure cuptions 

Fig.l. Giant compression and extension magneto-strain effects induced by magnetic field in Ni4 8 Mn 30 Ga22 alloy. 
Fig. 2. Different stages of twin boundary motion in situ observed in the rotating magnetic field 0.6T in Ni4 8 Mn 3 oGa22 
alloy. 

Fig. 3. Magnetization curves along easy and hard axes of Ni4 8 Mn3 Ga 2 2 sample constrained in a single variant state. 
Fig. 4. Magnetic anisotropy of Ni 48 Mn 3 oGa22 and field dependence of magnetization free energies for easy and hard 
magnetization directions 

Fig. 5. Two-variant twin microstructure, easy magnetization axes alignment and magnetic driving force responsible 
for twin boundary motion in a magnetic field applied. 

Fig. 6. Data of mechanical testing for Ni4 8 Mn3oGa22 sample during its transformation between two twin variants 
of martensitic phase. 

Fig. 7. Strain-stress behavior caused by transformation between two single twin variants induced by mechanical 
stress in Ni-Mn-Ga martensitc. 

Fig. 8. Model calculation result of the magnetic field induced strain effect in comparison with experiment. 

Fig. 9 Experimental data and model calculation of magnetization hysteresis caused by magnetic field controlled 
motion of twin boundaries in Ni-Mn-Ga. 
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